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ABSTRACT 

Context. The presence of pulsations in late-type Be stars is still a matter of controversy. It constitutes an important issue to establish 
the relationship between non-radial pulsations and the mass-loss mechanism in Be stars. 

Aims. To contribute to this discussion, we analyse the photometric time series of the B8IVe star HD 50209 observed by the CoRoT 
mission in the seismology field. 

Methods. We use standard Fourier techniques and linear and non-linear least squares fitting methods to analyse the CoRoT light 
curve. In addition, we applied detailed modelling of high-resolution spectra to obtain the fundamental physical parameters of the star. 
Results. We have found four frequencies which correspond to gravity modes with azimuthal order m = 0,-1, -2, -3 with the same 
pulsational frequency in the co-rotating frame. We also found a rotational period with a frequency of 0.679 cd~' (7.754 /jHz). 
Conclusions. HD 50 209 is a pulsating Be star as expected from its position in the HR diagram, close to the SPB instability strip. 

Key words. Stars: emission line, Be - Stars: oscillations (including pulsations) - Stars: individual HD 50 209 

1 . Introduction the same region of the HR diagram as f3 Cephei and SPB stars, 

it is generally assumed that pulsations have the same origin, i.e. 
Classical Be stars are defined as main sequence or slightly and/or g . mode pu i sat i ons driven by the K-mechanism associ- 

, evolved B-type stars whose spectrum has or had at some time ated wkh the Fe bump The current theoretical models have dif . 

g . one or more Balmer lines in emission. They are physically un- ficuhies in describing the pulsational characteristics of Be stars, 

derstood as rapidly rotating B-type stars with line emission aris- due tQ the high rotational velocity of these objects A brief de _ 

ing from a circumstellar disk in the equatorial plane, composed scription of ±e current kn ow i edge of th e pulsational b ehavior of 

of matter ejected from the stellar photosphere by mechanisms Be stars can be found in | Neiner eta i] rf20Q9h and lErruho et all 

not yet understood (see IPorter & Riviniusll2003L for a complete (2009) 

review). A significant fraction of Be stars show short-term pho- HD 50 209 is a late Be star of spectral t ype B8IVe and mag- 

tometnc and spectroscopic periodic variability, which is com- nitude y = g M ft hafj been smdied by Gutierrez . Soto et al . 

monly attributed to non-radial pulsations. As Be stars occupy ,(2007), using data from Hipparcos, ASAS-3 and the OSN 

(Observatorio de Sierra Nevada). The Hipparcos data analysis 

* Based on observations made with the CoRoT satellite, with FEROS yielded a f reqU ency of variation at 1 .689 cd 1 considered as un- 
at the 2.2m telescope of the La ^Sdla Observatory under the ESO Large and fr at L47 cd -i although with a 

Programme LP178.D-0361 and with Narval at the Telescope Bernard , , „ A^^t , i • i i Z 

t 7 f a tv a ^k a ■ r^u . lower amplitude. From OSN, the analysis revealed a frequency 
Lyot of the Pic du Midi Observatory. „n™ .-l r™ . . . , . „ J . „ _ , , , . 

HTU m . , i ui • i , f at 1.4889 cd . The analysis of the ASAS-3 dataset showed sig- 

Table |2] is only available in electronic form via J .„__ ,_i , , ,_i 

http://www.edpsciencies.org mficant peaks at frequency 2.4803 cd and 1 .4747 cd 



*** 



Current address: Valencian International University (VIU), Jose The occurrence of pulsations in late B-type stars has been a 

Pradas Gallen s/n, 12006 Castellon, Spain. matter of controversy in the recent literature. Hubert & Floquetl 

1 Current address: Laboratoire AIM, CEA/DSM-CNRS-Universite CHI) showed that pulsations in B6-B9 type s tars are much 

Paris Diderot; CEA, IRFU, SAp, centre de Saclay, F-91 191, Gif-sur- less common than in their early-type counterparts. Baade ( 1989) 

Yvette, France. failed to detect line profile variations in the spectra of B8-B9.5 
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stars. However. ISaio et aTl 12007) presented the detection of low 
amplitude g-modes in the B8 Ve star B CMi. To ascertain whether 
Be stars of all types do present pulsations is a key issue in order 
to establish the relation between non-radial pulsations and the 

mass ejection mechanisms. 

The CoRoT satellite (lAuvergne et all l2009t) observed HD 
50209 in its seismology field. The observations span 136 days 
in the Galactic anti-centre direction (LRA1), between October 
18th 2007 and March 3rd 2008, with a sampling of 32 seconds. 
The light curve contains 328 279 data-points with a duty cycle 
of 89%. 



2. Frequency analysis 

For the frequency analysis, we employed the code pasper 
(Dia go et al.l 120081), which is based on the cla ssical discrete 
Fourier transform dDeemingll 1 975b ISca rgle 1982) and the least- 
square fitting of a sinusoidal function in the time domain. When 
a frequency is found, it is prewhitened from the original data 
and a new step starts looking for a new frequency in the residu- 
als. The method is iterative and stops when the new frequency is 
not statistically significant. 

The criterion used to determine whether the frequencies are 
statistically significant is the signal to noi se amplitude ratio re- 
quirement described in Breg er et"ail(ll993l) . It consists of the cal- 
culation of the SNR of the frequency in the periodogram. This 
is made calculating the signal as the amplitude of the peak for 
the frequency obtained and the noise as the average amplitude 
in the residual perio dogram after the p rewhitening of all the fre- 
quencies detected. Breg er et"ail JT993) proposed that a value of 
SNR > 4 is a reliable criterion to distinguish between peaks due 
to real frequencies and noise. 

All the techniques used in this paper are described in detail 
by iGutierrez-Soto e t al. (2009). As mentioned in that paper, the 
orbital characteristics of the CoRoT satellite produces signal in 
the data at the orbital frequency (13.97 cd -1 , 161.689 //Hz) and 
day/night peaks (2.007 cd , 23.229 //Hz) and their harmonics. 
All those frequencies of instrumental origin have been removed 
from the list of detected frequencies. 

The re solution in frequency of our analysis is 1/2 T ~ 3.66 x 
10~ 3 cd" 1 ( Kalling^etaLl[2008 ), T being the total interval cov- 
ered by observations. The uncertainty on the detected frequen- 
cies is V3 x 10~ 6 cd~'. T his value has been derived analyticall y 
using the formula given by Montgomery & 0'Donoghue(1999), 
and takin g into account the correlations in the residuals, as de- 
scribed by ISchwarzenberg-Czernvl (1 1 99 ll) . 

In Fig. [I] (top panel) we present the CoRoT light curve of 
HD 50209. A long-term decreasing pattern is apparent. It is a 
common feature of the CoRoT light curves to show a linear or 
an almost linear decreasing pattern of instrumental origin due 
to the CCD ageing. On the other hand, most Be stars present 
long-term variations which in our case contribute to the observed 
variability. We have removed these trends by fitting a polynomial 
function of second degree without making assumptions on their 
origin, obtaining the light curve depicted in the bottom panel of 

Fig.LT] 

In Fig. [2] we display the Fourier spectrum, showing the 
typical aliasing at the o rbital frequency (as described in 
Gutierrez-Soto et al.l 120091) . From the frequency analysis, we 
obtain 60 significant frequencies which are listed in Table [2] 
(full table only available online). Note that frequencies around 
3/T ~ 0.022 cd _1 (0.255 //Hz) or lower correspond to periods 
of the order of the entire dataset time-span, and hence should 



Table 1. Principal frequencies of each set detected in HD 50 209. 



Set 


Freq. 
[cd- 1 ] 


Freq. 
[//Hz] 


Amp. 
[mmag] 


Main freq. 


Comment. 


Ft 


1.48444 
1.49028 
1.47860 
1.49613 


17.181 
17.248 
17.113 
17.316 


2.529 
1.616 
1.090 
0.653 


h 


/5+2/4 


F 2 


2.16238 
2.16822 


25.027 
25.095 


0.904 
0.495 


h 


fs +3/ 4 


F 3 


0.79482 
0.77874 
0.80650 


9.199 
9.013 
9.334 


0.607 
0.407 
0.358 


h 


fs + fi 


Fa 


0.67939 
0.69108 


7.863 
7.998 


0.425 
0.368 


U 


fro, 


Fs 


0.10811 


1.251 


0.202 


h 




F 6 


2.96889 


34.362 


0.121 


h 
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Fig. 1. Top: CoRoT light curve of the star HD 50209. Bottom: 
CoRoT light curve detrended with a polynomial of second de- 
gree. 
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Fig. 3. The 60 significant frequencies detected for the star HD 
50 209 revealing six different and separated frequency groups. 
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Fig. 2. Left: Complete Fourier spectrum for HD 50209 showing aliasing effects at the orbital frequency 13.97 cd . Right: Fourier 
spectrum for HD 50 209 in the domain of the detected stellar frequencies. 



Table 2. Complete list of significant frequencies obtained with pasper for the star HD 50209. The time reference for the phase is 
HJD = 2454391.95. The full table is available online. 



# 


Freq. 


Freq. 


Amp. 


Amp. error 


Phase 


Phase error 


S/N 


Comment. 




[cd" 1 ] 


L"Hz] 


[mmag] 


[mmag] 


[0,1] 


[0,1] 






1 


1.48444 


17.181 


2.5299 


0.0011 


0.7348 


0.0001 


265 


fi 


2 


1.49028 


17.2486 


1.6167 


0.0012 


0.1577 


0.0001 


169 




3 


1.47860 


17.1134 


1.0901 


0.0011 


0.2747 


0.0002 


114 




4 


2.16238 


25.0275 


0.9043 


0.0010 


0.9746 


0.0002 


94 


h 


5 


0.79482 


9.19931 


0.6077 


0.0011 


0.9866 


0.0003 


63 


h 



be considered with caution. They might be produced by the de- 
trending process or artifacts of the frequency analysis. However, 
we have included them for completeness. 

The entire set of detected frequencies together with their 
amplitudes is plotted in Fig. [3] The frequencies found are dis- 
tributed in 6 main groups (F, with i = 1, ...,6). The frequencies 
with the highest amplitudes are listed in Table [TJ All groups are 
clearly separated, except those centered on fa and ft,. In each 
group we note the existence of equidistant multiplets separated 
by intervals marginally larger than our frequency resolution. In 
Fig.|4]we provide the diagrams of the folded light curve with the 
two main frequencies f\ and fi- 

In addition, we performed a time-frequency analysis of the 
original light curve by applying the pasper code to sliding 
windows with different sizes (for details on the method, see 
iHuat e t al. 2009). The results presented in Fig.|5]show that all the 
main frequencies show large amplitude changes. See Section [4] 
for a detailed discussion of these changes. 



3. Ground-based observations 

High resolution spectroscopy of HD 50 209 was obtained with 
the FEROS spectrograph at the 2.2m telescope in La Silla, 
in the framework of a large program complementary to the 
CoRoT mission. Observations were carried out one year prior 
to the LRA1 run of CoRoT to investigate spectroscopically 
the rapid variability prev iously detected in photometry by 
Gutierrez-Sot o et ail (120071) . Seventy spectra with high spectral 
resolution and high signal-to-noise ratio were obtained. Fourteen 
additional spectra were obtained with the Narval spectropo- 
larimeter at the TBL at Pic du Midi Observatory (France), con- 
temporary to the CoRoT run. However, due to bad weather con- 
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Fig. 4. Top: Phase diagram for the HD 50 209 data folded with 
the frequency ft = 1.48444 cd _1 (17.181 //Hz). Bottom: Phase 
diagram for the HD 50 209 data, prewhitened from the frequency 



ft and folded with the frequency ft = 2.16238 cd 
//Hz). Note that two cycles are shown for clarity. 



(25.027 



ditions, the signal-to-noise ratio is not high enough to search for 
low amplitude variability in the Narval data. 

FEROS observations (R = 48000, 3600 - 9200 A) were re- 
duced with MIDAS (wavelength calibration, bias and flat field 
corrections and earth motion correction). For the Narval obser- 
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Fig. 5. 3-D view of the frequencies determined using the time- 
frequency analysis with a sliding window of 30 days. The color 
scale indicates the amplitude of the detected frequencies in HD 
50209. The time origin corresponds to HJD = 2454391.95. 



vations (R = 65000) the sum of the 4 sub-exposures obtained in 
Stokes V sequences were used. Exposures w ere reduced locally 
with LibrEsprit based on the Esprit software dDonati et al.f l997) 
and then summed. The continuum normalization was carried out 
with IRAF. 

The fundamental physical parameters of HD 50209 have 
been accurately determined from the newly avail able spectro- 
scopic data adopting a procedure described by Frema fet al.l 
2006) and corre cting for gravitational darkening effects 



(Fre mat et alj |2005). The obtained parameters are derived by fit- 
ting to the spectra, adopting the different values of Q/£2 e pre- 
sented in the first row of Table [3] The results are presented in 
Tablela nd are consisten t with the previous parameter determi- 
nation dFremat et al.ll2006l) . 

Spectroscopic data also show that emission is present in the 
first Balmer lines and visible up to H6. Ha is a strong emission 
line (EW = —19 A , Imax = 4.68 Icont) with inflection points 
on the wings and close double peaks at the centre. The line pro- 
file i s typical of a Be star s een under a moderate inclination an- 
gle dHanuschi k et al.|[l996l) . HB, Hy and H5 show a symmetrical 
and double-peaked emission embedded in the broad wings of the 
photospheric profile. Moreover, numerous Fe n and Cr n lines as 
well as the IR Ca n triplet and O i 8446 A lines are in emission 
with a double component. Note that we have not found signifi- 
cant changes in the global emission state of HD 50 209 between 
the FEROS and Narval spectra taken about one year apart. 

In addition, forbidden single-peaked [O i] 6300 A and [Fe n] 
emission lines are detected with a very weak intensity (/ < 
1.01% of the continuum, see Fig. |6|. The star is isolated, out- 
side any formation region, and without the large IR excess typ- 
ic al of HAeBe stars. I t coul d be an extreme Be star according 
to lde Winter & Pere3 d 19981) or a n unclB[e] star in the scheme 
proposed by Lame rs et al.l dl998), though the B[e] character is 
very weak. 

The mean Ha, HB and Fe II 5 197.6 A and 53 16.6 A profiles 
are shown in Fig. [7] The separation between the emission peaks 
and the lack of a shell profile exclude an inclination angle close 
to 90 degrees, and hence the physical parameters in the first row 
of Table [3] are not suitable. This implies that the rotational fre- 
quency is at least 90% of the critical velocity. The shape of the 
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50209. 



5 
4 
3 
2 

'55 1 
I 



.51.1 

DC 1 

0.9 

0.8 
0.7 



Ha 
















6555 


6560 


6565 


6570 




4855 


4860 


4865 


4870 



1.06 
1.04 
1.02 
1 

0.98 

1.08 
1.06 
1.04 
1.02 
1 

0.98 





f\ Fell 5197.6 ' 


5190 5195 


5200 5205 


Fell 5316.6 / 





5310 5315 5320 5325 



Wavelength [A] 

Fig. 7. Left panels: Ha and HB emission. Right panels: Different 
emission lines in the Fen region in the star HD 50209. In all 
panels the horizontal axis is wavelength and the vertical one is 
relative intensity. 



Ha line profile and the forbidden lines discussed above indicate 
the presence of a very extended circumstellar disk. 

In spite of the high quality of the FEROS spectra and the 
presence of rapid photometric variability apparent in the CoRoT 
data (see Fig. [TJ, we have not been able to detect periodic vari- 
ability in purely photospheric lines of HD 50 209. Some variabil- 
ity is present slightly above the detection limit but no consistent 
periodic behavior could be established. Weak rapid variations 
were detected in the intensity of the V and R emission com- 
ponents of Balmer lines and their V/R ratio. Note also that the 
quantity (V + R)/2 seems to be slightly lower in the second part 
of the ESO run, suggesting a slight weakening of the emissivity 
of the disk. As shown in the mean variance in the HB profile, 
the variability is concentrated in the emission part. Note that the 
variability is more conspicuous on certain days. Unfortunately, 
the gap between the two parts of the ESO run, the alternating 
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Table 3. Fundamental parameters for HD 50209 computed for different values of Cl/Q. c and corrected for veiling. The errors in 
T e ff, logg and V sin i are taken from the apparent parameters. 



n/n c 


i 




Tss 




l°gg 


V sin i 


M 


L 






frol 




[dej 


il 


[K] 




[cgs] 


[km/s] 


[M e ] 


[io] 


[Ra] 


[cd~ 


] 


0.80 


89 ± 


20 


13600 ± 


1500 


3.60 ±0.11 


192 ± 20 


4.67 ± 0.95 


2.97 ± 


0.39 


6.44 ± 1.53 


0.60 ± 


0.20 


0.90 


64 ± 


14 


13600 ± 


1500 


3.56 ±0.11 


205 ± 20 


4.76 ± 1.02 


3.02 ± 


0.39 


7.27 ± 1.74 


0.62 ± 


0.20 


0.95 


57 ± 


15 


13600 ± 


1500 


3.53 ±0.11 


211 ±20 


5.07 ± 1.15 


3.08 ± 


0.40 


8.06 ± 2.03 


0.62 ± 


0.22 


0.99 


55 ± 


12 


13400 ± 


1500 


3.46 ±0.11 


221 ± 20 


4.99 ± 1.22 


3.14 ± 


0.41 


9.34 ±2.35 


0.58 ± 


0.20 



Frequency [uHz] 

16.666 17.129 17.593 

' ' Periodogram for f 1 

Spectral window for HD 50209 




1.42 1.44 1.46 1.48 1.5 1.52 1.54 

Frequency [c/d] 

Fig. 8. Periodogram of the original light curve (solid line) and 
spectral window shifted at fa = 1.48444 cd -1 (dashed line). Note 
that the full width at half maximum of the peak in the peri- 
odogram at the frequency fa is broader than expected from the 
spectral window, suggesting the presence of more frequencies or 
amplitude changes. The vertical line indicates the position of the 
frequency fa . 



of day/night and the weakness of variations prevent any reliable 
estimate on the time-scale of variations. 



4. Discussion 

As a result of our analysis, we have obtained 60 significant fre- 
quencies, grouped in six sets. In each set, different equidistant 
frequencies close to the main detected ones are present. The 
Fourier analysis of the light curve by means of sliding windows 
results in the detection of a significant variation of amplitude of 
the main frequencies (Pig. [5j. 

The first issue to be addressed is the presence of frequency 
multiplets and amplitude variations. Both phenomena are re- 
lated, and can be due either to the presence of close frequencie s 
or actual amplitude variations dBreger & Pamvatnvkhl |2006). 
The beating of two or more close frequencies will appear as 
a single frequency with variable amplitude when we split the 
data sample in shorter intervals, with the consequent loss of fre- 
quency resolution. On the other hand, true amplitude variability 
will produce peaks in the power spectrum broader than expected 
from the frequency resolution (Fig. [8), and the prewhitening of 
the main frequency will leave power in the wings of the main 
peak which will lead to false frequency doublets or multiplets. 

In order to discriminate between the beating of several fre- 
quencies and an amplitude chan ge of a single frequency we ha ve 
applied the method described in lBreger & Pamvatnvkhl d2006l) to 
study the relationship between amplitude variability and phase 



variations of an assumed single frequency. The idea behind this 
method is that the beating of frequencies will produce a phase 
variation with the same period as the beating period, while in 
the case of a single frequency with variable amplitude, the phase 
will remain constant. 

We studied the six frequency groups with the abovemen- 
tioned method. However, the results obtained are inconclusive, 
for two main reasons: i) the beat periods of the close frequencies 
are larger than the time coverage of CoRoT data, and hence, we 
cannot evaluate the consistency of the phase variations; ii) the 
phase variations predicted by the models with four or more fre- 
quencies are very small. At our detection level we were not able 
to discriminate between the predicted low amplitude variations 
and no variation at all. Consequently, we cannot firmly reject 
either of the two interpretations. 

In the following, we analyze the six main frequencies found, 
which, as discussed above, can be either single frequencies with 
variable amplitude or groups of close frequencies around a cen- 
tral value. The highest frequency, fa, is the first harmonic of fa, 
i.e. fa =2f\, and hence it will not be considered in the analysis. 
Frequencies fa, fa, fa and fa are equidistant within the frequency 
resolution at the 3-cr level, and the separation between them is 
the frequency fa. Moreover, fa = 0.61939 cd~'is consistent with 
the rotational frequency of the star as given in Table [3] As a con- 
sequence, we have only two independent main frequencies, fa 
and fa, and all the others are of the form /; = fa + nfa with 
n= 1,2,3. 

Let us recall that the observed frequencies in a rotating star 
are related to the pulsational frequency in the co-rotating frame 
by the expression: 

v= \v-mQ\ (1) 

where v is the frequency in the co-rotating frame and D. is the 
rotational frequency of the star. From this expression, and con- 
sidering fa as the rotational frequency as discussed above, the 
four frequencies fa, fa, fa and fa can be consistently interpreted 
as modes with the same pulsational frequency in the co-rotating 
frame and with m = 0,-1, -2, -3 respectively. 

The fundamental parameters presented in Table Opiace HD 
50 209 marginally o utside the SPB instability strip calculated by 
iPamyatnykhl (LL999), although its position is compatible with the 
strip at the 1-cr error level. Hence, we consider that the pulsations 
detected are gravity modes typical of S PB stars. HD 50 209 is a 
SPBe star, the designation proposed bv lWalker etal.1 d2005l) for 
the Be stars pulsating in g-modes, considered as rapidly rotating 
counterparts of the SPB stars. 

Due to the fact that in late-type B stars the frequencies of 
the g-modes in the co-rotating frame are much smaller than 
the rotational frequency, the frequencies of these modes in the 
observer's frame are close to |raQ|. This leads to the diffi- 
culty that the observed frequencies close to the expected rota- 
tional frequency ca n be either interpreted as g- mode pulsation s 
dWalker et alj |2005) or as rotational modulation (Bal onalll995l) . 
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In our case, we can for the first time discriminate between the 
rotational frequency and the frequency of the g-mode pulsation 
with azimuthal order \ m\ = 1 , as we have detected both of them 
clearly separated by an interval much larger than the frequency 
resolution, thanks to the long duration of observations and pre- 
cise measurements of CoRoT. We can be confident that the fre- 
quencies /s, f}, f\ and ft are not related to the rotation and hence 
they are true gravity mode pulsations. 

The frequency in the co-rotating frame of the detected 
modes, namely f$ = 0.10811 cd~', is significantly lower than 
the frequencies commonly found for both SPB and SPBe stars. 
However, this low value is consistent with what is expected for 
an m — mode in models of rapidl y rotating late-type B stars 
(IWalker et alJl2005t ISaio et al.ll2007h . The modelling of the de- 
tected pulsations will provide more insights into the nature of 
the pulsational modes. 

5. Conclusions 

The high precision photometry and long duration of continuous 
observations provided by the CoRoT mission has allowed the de- 
tection of g-mode pulsations in the late-type Be star HD 50 209. 
This supports the fact that all Be stars have non-radial pulsations 
that could play a critical role in the mass ejection mechanism. 

From our analysis, we have found pulsation in four modes 
with the same frequency in the co-rotating frame with azimuthal 
order m = 0,-1,-2,-3. We have also detected the rotational 
frequency, both as a significant peak in the power spectrum and 
as the separation of frequencies with different m. The accurate 
determination of the rotational period will play an important part 
in constraining the fundamental parameters of the star in order 
to perform the seismic modelling. 

For the first time, we have been able to observe simultane- 
ously the rotational frequency and the pulsational frequencies 
and separate them, implying that the frequencies we attribute to 
g-mode pulsations cannot be interpreted as the effect of the ro- 
tational modulation. This constitutes a proof of the presence of 
pulsations in HD 50209. 
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Table 2. Complete list of significant frequencies obtained with pasper for the star HD 50209. The time reference for the phase is 
HJD =2454391.95. 



# 


Freq. 


Freq. 


Amp. 


Amp. error 


Phase 


Phase error 


S/N 




[cd- 1 ] 


L"Hz] 


[mmag] 


[mmag] 


[0,1] 


[0,1] 




1 


1.48444 


17.181 


2.5299 


0.0011 


0.7348 


0.0001 


265 


2 


1.49028 


17.2486 


1.6167 


0.0012 


0.1577 


0.0001 


169 


3 


1.47860 


17.1134 


1.0901 


0.0011 


0.2747 


0.0002 


114 


4 


2.16238 


25.0275 


0.9043 


0.0010 


0.9746 


0.0002 


94 


5 


0.79482 


9.19931 


0.6077 


0.0011 


0.9866 


0.0003 


63 


6 


1.49613 


17.3163 


0.6531 


0.0011 


0.2873 


0.0003 


68 


7 


2.16822 


25.0951 


0.4953 


0.0010 


0.2636 


0.0003 


51 


8 


0.01461 


0.16909 


0.4324 


0.0010 


0.2149 


0.0004 


45 


9 


0.77874 


9.01319 


0.4074 


0.0010 


0.0079 


0.0004 


42 


10 


0.80650 


9.33449 


0.3583 


0.0010 


0.4104 


0.0005 


37 


11 


0.67939 


7.86331 


0.4256 


0.0010 


0.3579 


0.0004 


44 


12 


0.69108 


7.99861 


0.3685 


0.0009 


0.5289 


0.0004 


38 


13 


0.82258 


9.5206 


0.2768 


0.0010 


0.4198 


0.0006 


29 


14 


0.80066 


9.2669 


0.2740 


0.0011 


0.1969 


0.0006 


28 


15 


0.10811 


1.25127 


0.2030 


0.0010 


0.5750 


0.0008 


21 


16 


0.81527 


9.436 


0.1725 


0.0009 


0.2895 


0.0009 


18 


17 


0.78897 


9.1316 


0.1214 


0.0011 


0.8024 


0.0014 


12 


18 


0.09350 


1.08218 


0.1509 


0.0009 


0.4966 


0.0010 


15 


19 


0.03506 


0.40578 


0.1359 


0.0009 


0.3221 


0.0011 


14 


20 


0.11396 


1.31898 


0.1647 


0.0010 


0.7770 


0.0010 


17 


21 


1.51074 


17.4854 


0.1217 


0.0009 


0.2095 


0.0012 


12 


22 


0.65894 


7.62662 


0.1333 


0.0010 


0.2271 


0.0011 


13 


23 


2.96889 


34.3622 


0.1210 


0.0009 


0.7144 


0.0012 


12 


24 


0.07451 


0.86238 


0.0862 


0.0009 


0.2504 


0.0017 


9 


25 


0.64725 


7.49132 


0.0980 


0.0009 


0.1021 


0.0015 


10 


26 


0.82842 


9.58819 


0.1421 


0.0010 


0.7036 


0.0011 


14 


27 


0.00876 


0.10138 


0.1040 


0.0010 


0.6659 


0.0015 


10 


28 


0.11980 


1.38657 


0.1164 


0.0010 


0.9496 


0.0013 


12 


29 


0.83573 


9.6728 


0.1160 


0.0010 


0.5225 


0.0013 


12 


30 


1.50197 


17.3839 


0.1158 


0.0011 


0.4101 


0.0014 


12 


31 


0.84157 


9.74039 


0.1131 


0.0010 


0.7923 


0.0014 


11 


32 


0.66916 


7.74491 


0.0847 


0.0009 


0.0830 


0.0018 


8 


33 


0.74952 


8.675 


0.0751 


0.0009 


0.0087 


0.0019 


7 


34 


0.04967 


0.57488 


0.1041 


0.0009 


0.5844 


0.0014 


10 


35 


0.70131 


8.11701 


0.0854 


0.0009 


0.8694 


0.0017 


8 


36 


0.02337 


0.27048 


0.0867 


0.0009 


0.2750 


0.0017 


9 


37 


0.06136 


0.71018 


0.0894 


0.0009 


0.8152 


0.0017 


9 


38 


0.56397 


6.52743 


0.0792 


0.0009 


0.1145 


0.0018 


8 


39 


1.45668 


16.8597 


0.0774 


0.0009 


0.4253 


0.0019 


8 


40 


0.62533 


7.23762 


0.0694 


0.0009 


0.1556 


0.0021 


7 


41 


0.54205 


6.27373 


0.0726 


0.0009 


0.8328 


0.0020 


7 


42 


0.86203 


9.9772 


0.0731 


0.0009 


0.7081 


0.0020 


7 


43 


2.15361 


24.926 


0.0749 


0.0009 


0.2862 


0.0019 


7 


44 


0.72030 


8.33681 


0.0829 


0.0009 


0.8066 


0.0018 


8 


45 


2.96158 


34.2775 


0.0672 


0.0009 


0.4195 


0.0021 


7 


46 


0.73199 


8.47211 


0.0681 


0.0009 


0.0589 


0.0022 


7 


47 


0.12711 


1.47118 


0.0683 


0.0009 


0.1174 


0.0021 


7 


48 


2.97619 


34.4466 


0.0640 


0.0009 


0.0250 


0.0022 


6 


49 


2.17406 


25.1627 


0.0875 


0.0010 


0.5229 


0.0019 


9 


50 


2.17991 


25.2304 


0.0755 


0.0010 


0.7423 


0.0021 


7 


51 


1.53996 


17.8236 


0.0590 


0.0009 


0.4461 


0.0024 


6 


52 


0.77290 


8.9456 


0.0615 


0.0010 


0.8969 


0.0025 


6 


53 


1.38363 


16.0142 


0.0553 


0.0009 


0.7286 


0.0026 


5 


54 


0.47484 


5.49583 


0.0574 


0.0009 


0.4781 


0.0025 


6 


55 


0.61510 


7.11921 


0.0560 


0.0009 


0.2581 


0.0026 


5 


56 


1.52535 


17.6545 


0.0559 


0.0009 


0.2121 


0.0026 


5 


57 


0.87079 


10.0786 


0.0642 


0.0009 


0.8295 


0.0023 


6 


58 


0.88394 


10.2308 


0.0567 


0.0009 


0.9051 


0.0026 


5 


59 


0.58150 


6.73032 


0.0525 


0.0009 


0.6858 


0.0028 


5 


60 


0.17971 


2.07998 


0.0515 


0.0009 


0.1046 


0.0028 


5 



S/N Comment. 
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